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We consider a model in which Goldstone bosons, described by a SU(N) chiral nonlinear a model, 
are coupled to an 7V-plet of colored fermions by means of a Yukawa interaction. We study the 
one-loop renormalization group flow and show that the non-Gaussian UV fixed point, which is 
present in the purely bosonic model, is lost because of fermion loop effects unless N is sufficiently 
large. We then add four-fermion contact interactions to the lagrangian and show that in this case 
there exist several non-Gaussian fixed points. The strength of the contact interactions, predicted 
by the requirement that the theory flows towards a fixed point in the UV, is compared to the 
current experimental bounds. This toy model could provide an important building block of an 
asymptotically safe model of the weak interactions. 
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I. MOTIVATIONS 

o 

Although the nonlinear a model (NLcrM) is usually seen — because of its perturbative nonrenormalizability — as 
a mere low-energy effective field theory, it could also be renormalizable in a nonperturbative sense, namely at a 
nontrivial (non-Gaussian) fixed point (FP) of the renormalization group (RG) flow [1] (see [2] for the application of 
this idea in the context of gravity); in this case, the model is said to be asymptotically safe (AS). 

As appealing as this approach may seem, it has the disadvantage that perturbation theory can only be a rough 
guide, and therefore many results are necessarily only suggestive rather than conclusive. If in spite of this serious 
drawback one considers the RG one-loop /3-functions, or some resummation thereof, it is easy to see that a nontrivial 
r— i FP is indeed present in the NLcrM [31 H] - It persists when one couples the Goldstone bosons to gauge fields [5] and, in 
some cases, when one adds terms with four derivatives of the Goldstone bosons [B]. While these results do not prove 
the existence of the FP, they are suggestive and we use them to justify our rather phenomenological approach in which 
O ! the one-loop results are assumed to hold at least qualitatively in the non-perturbative solution and the consequences 
1} of such an assumption are then worked out. 

Whether the physical system described by the NLcrM is actually on a RG trajectory leading to such a FP or not 
depends on the initial conditions at some given energy scale and the experimental constraints on the coefficients in 
| the operator expansion. Whereas it appears that this is not the case for the chiral NLcrM of strong interactions (for 
which the Goldstone bosons are the pions), this is still an open question in the case of the Goldstone bosons of the 
spontaneous breaking of electroweak interactions. This question will eventually be answered by the experimental 
data. 

In this work we address a more specific problem and reconsider the chiral NLcrM to study what happens to the FP 
of Rcf. [4 when the Goldstone bosons are coupled to fermions. This is an important question in view of utilizing the 
model in a realistic theory of the weak interactions. 

As we shall see, a one-loop computation shows that the inclusion of fermion fields drastically modifies the asymptotic 
properties of the NLcrM. The Yukawa coupling does not have a UV fixed point and instead grows with the energy. 
The Goldstone boson coupling hovers about the FP of the uncoupled model — as long as the Yukawa coupling size is 
small enough to be negligible — but is eventually pushed away by the increasingly large Yukawa coupling. 

We interpret this instability as a signal of new physics associated to the standard model (SM) fermions, and in 
particular the top quark, and model it by means of four-fermion contact interactions — somewhat along the lines 
of what is done in nuclear physics in describing the nucleon potential. The fermionic content of the model is left 
unchanged and identical to that of the SM. The model thus extended admits a number of non-trivial FPs and can 
thus be used in modeling the interaction between Goldstone bosons and fermions in a AS theory. 

We also show how current experimental bounds on quark contact interactions already limit the choice of RG 
trajectories which can be used. As the bounds are improved, or the new interactions are discovered, they will provide 
a direct test of the model. 
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II. FERMIONS AND GOLDSTONE BOSONS 



We consider a ST/(iV)-valued scalar field U — exp(j/7r a T a ) where ir a are the Goldstone boson fields, trT a Ti> = 5 a b/2 
and / is the Goldstone boson coupling, which in the SM case can be identified with 2/v, where v = 246 GeV is the 
Higgs VEV. We limit ourselves to the lowest order term in the NLctM lagrangian, which reads 

C a = ~Tr (U^d^UU^d^U) . (1) 

This model is invariant under separate SU(N)l and SU(N)r transformations acting on U by left- and right multi- 
plication respectively. 

We couple the Goldstone bosons to left- and right-handed fermions and ip^ carrying the fundamental represen- 
tation of SU(N)l and SU(N)r respectively (corresponding to the indices i = 1, . . . ,N), and also the fundamental 
representation of a color group SU(N C ) (corresponding to the indices a — 1, . . . ,N C ). In the real world the latter 
group is gauged; here we merely retain it as a global symmetry to count fermionic states. We couple the fermions in 
a chiral invariant way to the U field by adding to the NLerM lagrangian the fermion kinetic and the Yukawa terms: 

C. r = fcvfd^L + fai-fd^R - j ffiU**^ + h.c.) , (2) 

where we have explicitly written out the group indices in the interaction. 

In a more realistic model, the Yukawa coupling should distinguish among the fermion components and there should 
be more than one family, as it is the case in the SM, but this is not essential for the present discussion. Indeed, we find 
that the only phenomenologically relevant fermionic contribution comes from the top quark, so h in ^ will be viewed 
as the top Yukawa coupling. We also neglect the strong, weak and electromagnetic gauge couplings (after having 
checked that they do not alter our conclusions) and for this reason the derivatives in eqs. ((T])-(|2| are not covariant. 
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FIG. 1: Running of / and h (rescaled by a factor of 10 to fit the figure) for N — 2, N c — 3. The AS behaviour of / (represented 
in the figure by the blue dashed line) is destabilized around t — 4.5 (that is, around 22 TeV) due to the increasingly large 
Yukawa coupling contribution (green dashed line) which at about the same value becomes strongly coupled, that is larger 
than 2tt. The asymptotically safe behavior (continuous lines) is recovered after the introduction of the four-fermion contact 
interactions, as discussed in section III. 



Using a sharp cutoff regularization, we obtain the following one-loop RG equations for / — kf and h: 

— = * 4AT c -2 /i 3 H hf (4) 

dt 16tt 2 \ N ) 64tt 2 N J 1 ' 

The variable t is defined as logfc/fc where k is the renormalization group energy scale and k is a conventional 
reference energy, which we take to be ko = v. 
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The system of equations admits a number of possible UV fixed points. The Gaussian FP is the solution 

/ = 0, h = 0. This FP has one irrelevant (IR attractive) direction / and one marginal direction h. Any RG trajectory 
starting from it in the UV would be trivial at any scale. We therefore reject this choice because it is uninteresting. We 
emphasize that conventional chiral perturbation theory is the study of small deformations of the FP f — 0, interpreted 
as an IR FP. 

There is also a nontrivial FP at = 0, /* = Sw/yN for which / is a relevant (UV attractive) direction and h is 
marginally irrelevant. Again, interpreting such a FP as an UV attractive FP would imply the triviality of the Yukawa 
coupling at all scales. We thus see that this FP is the same as in the purely bosonic theory [4]. 

A physically viable UV limit for ([3])- ([4]) requires the existence of aFP with both h and / being nonvanishing. If 
h is treated as a i-independent constant, the /3-function for / has a zero at /* = 4^(471-2 +N c h 2 )/N, which is a 
deformation of the one appearing in the pure bosonic model. The existence of a nontrivial FP for the coupled system 
thus hinges on the existence of a nontrivial zero in the /3-function of h. This requires that the first term in the right 
hand side of eq. Q be negative, which is true for N > 2N C . 

Unfortunately this condition is not satisfied for the phenomenologically most important case N = 2, N c = 3. This 
is illustrated by the dashed curves in Fig. [l] for the initial condition f = 2 and ho = 0.7 at t = (thus mimicking 
the top-quark value). The first term on the right hand side of eq. ^ is initially dominant, leading to linear growth 
of /. The second term then grows in absolute value and at some point nearly balances the first one, leading to an 
approximate FP behavior in some range of energies. Eventually h, whose /3-function is everywhere positive, becomes 
large and the third term dominates leading to a Landau pole and the loss of AS. The scale at which destabilization 
occurs is very sensitive to the initial conditions and for the Yukawa couplings corresponding to light fermions no 
destabilization takes place up to very large energies. 

We conclude that our model is not AS in the case N — 2, N c — 3, in the one loop approximation. For it to be AS, 
either the one loop approximation must break down, or else new physical effects must enter in the fermion sector at 
some energy scale. 

It is interesting to compare this behavior to similar models. 

If the color symmetry was gauged there would be an additional contribution to eq. Q: 

The system ([3|-(|4]) would then resemble the SM where the second term in eq. Q is abstent. This model is found to 
predict a perturbative Yukawa coupling up to very high scales, in agreement with [7]. On the other hand, in the case 
of the NLcrM the second term in eq. Q is present and the qualitative behavior shown in Fig. [I] would not be changed 
by the term (J5J. 

A strictly related model is the linear sigma model coupled to one right-handed and Nl left-handed fermions, 
studied in (jjjjj. Our Goldstone modes are contained in their scalar sector, with the VEV v = 2// corresponding to the 
minimum of the scalar potential. There it was found that the the scalar potential and the Yukawa coupling admit a 
FP for 1 < Nl < 57. Our results are quite different due both to the different fermion content and to the non-linear 
boson-fermion coupling. 

III. FOUR-FERMION INTERACTIONS 

New physics associated to the SM fermions might restore asymptotic safety. In this section we will discuss a class 
of nontrivial UV FPs with asymptotically free Yukawa couplings that emerge once short-range interactions among 
the fermions are included. 

We now restrict ourselves to the case N = 2 and add to the lagrangian a complete set of SU (2) l x SU (2) r invariant 
four fermion operators. Requiring P invariance, all possible chiral invariant operators, up to Fierz reorderings are 
given by the following lagrangian: 

+ A 3 (^7^MV4 6 +^7^fl^V^?) 

+ A 4 (ffi^%fifrK+$R^Rp^1%) ■ (6) 

The coefficients have inverse square mass dimension. We do not include in the lagrangian in eq. ^ operators 
defined by taking the square of the Yukawa term in eq. ([2| because they are higher order from the point of view of 
chiral perturbation theory. 
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Strictly speaking, only the third SM fermion generation requires new physics to emerge at relatively low scales, since 
four-fermion operators involving the first two generations can be suppressed by much larger scales without spoiling 
our AS scenario. Yet, in discussing the experimental bounds on our model, see Section IV, we will be as conservative 
as possible and consider the rather pessimistic scenario in which |6]) also involve the first generation. We thus tacitly 
assume that the operators in Q are consistent with FCNC bounds. 

The symmetries imposed on eq. ^ make this lagrangian the minimal choice, and the one we will study here. 
More general sets of operators may well be relevant depending on the symmetries of this new sector; the SM group 
SU(2) x U(l) being the first instance coming to mind. The analysis is in this case more cumbersome because more 
operators must be included, but no distinctive features are expected to arise. 

The operators in eq. ^ are similar to those discussed in top-quark condensation models. In these and other models 
of composite quarks only operators with vector current structure that is iso- and color singlet are usually considered. 
Here the full set of operators are taken into account since their couplings mix in the RG evolution equations. Notice 
however that the RG regime we are interested in is very different because we do not seek to model chiral symmetry 
breaking and therefore the IR initial conditions for the parameter Ai are taken at finite values and no phase transition 
is present. A discussion of the four-fermion lagrangian in eq. ([6]), and its FP, as a model of chiral symmetry breaking 
can be found in [3]. The special role played in that context by the operator associated to Ai has been recently 
emphasized in 10J . 

We write the system of coupled /3-function equations for the dimensionless variables /, h and Ai = \k 2 . In these 
variables, we find the following RG equations: 
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In the RG equations above for the coefficients A.; we have neglected contributions coming from the Yukawa terms 
which are proportional to h 2 f 2 , h 2 Xi or h 2 X 2 /f 2 . These latter terms are negligible in the UV because we will select 
fixed points for which the Yukawa coupling approaches zero. They are also subdominant in the IR because the 
couplings A 2 , / 2 are IR-free, and such that A 2 // 2 — > at low scales. 



Notice that only the operators proportional to Ai and A 2 contribute to the /3-function for h. The other two operators 
do not contribute because of their chiral properties. Moreover, it is crucial that the combination 16(-/V c Ai + A 2 ) at 
the FP of Ai and A 2 be different from zero and larger than f 2 because otherwise there would be no UV FP. 

From here on we fix JV C = 3. We begin by observing that the /3-functions of the Ai form a closed sub-system. 
The numerical study of these equations reveals the presence of 16 real fixed points. Their coordinates are given in 

the first four columns in Table I. They are listed in order of decreasing trace of the stability matrix -t^-S from the 

most UV-repulsive, the Gaussian FP fpO, to the most UV attractive fp4. The coefficients of the operators related 
to UV-repulsive directions are completely determined by the AS condition. Hence, the FPs with a low number of 
UV-attractive directions are the most predictive ones, and hence the most phenomenologically appealing. The number 
in the name of the FP is the number of relevant (UV attractive) directions. These values can then be used to find 
the zeroes of f3h, and these in turn are used to find the zeroes of ftj-. 

We are interested in the sixteen FPs for the complete system with = 0; this requirement implies f 2 — 32tt 2 . At 
each of these sixteen FPs the direction / is always a relevant one, with eigenvalue —0.45; the direction h is also an 
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0.5 


fpla 





-28.71 


-7.18 





1.22 


fplb 








7.85 


-9.51 
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TABLE I: Values of the coefficients A« for the 16 FPs discussed in the text. /» = 17.78, h* = for all FPs. The FP fplc is 
boxed. 
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FIG. 2: Running of the Ai for the FP fplc. Ai and A4 are equal to zero at all energies. 



eigendirection, with eigenvalue 

(/?-16(JV c Ai. + A 2 .)) • (8) 

The numerical values of are listed in the last column of Table I. The seven FPs with eh > are physically 
uninteresting since the requirement of flowing to one of them in the UV implies that h(t) = at all scales. The other 
nine FPs of the fermionic sector for which < also admit a FP with nonzero h. We are not interested here in these 
additional nontrivial FPs, but restrict our attention to those with h* = 0. For the A^, the condition of flowing to fpnx 
in the UV yields 4 — n predictions. The values of / and h remain always free parameters, to be fixed by comparison 
with the experiment. 

We have studied numerically the trajectories emerging from the FPs in the directions of the relevant eigenvectors. 
Some of them lead to divergences, others flow to other FPs. Here we shall consider only the single renormalizablc 
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trajectory that ends at fplc in the UV. This is a natural choice as it is the most predictive (the one with the smallest 
number of relevant directions) among the FPs with < 0. 

In order to select the fine-tuned initial conditions in the IR that guarantee AS at fplc, we have first solved numerically 
the flow equations of the fermionic subsystem for decreasing t, starting at an initial point Aj* + lCP 8 «i, where Vi is the 
relevant eigenvector. This trajectory is attracted towards fpO after roughly 20 e-foldings, so the four-fermion couplings 
can be taken arbitrarily small by selecting the IR value of t appropriately. We then shift t such that this value is zero, 
in accordance with our convention that t = corresponds to the scale k — v. Now we pick a trajectory for the whole 
system by fixing the initial values of the Xi to agree with the ones find by this method, while the initial value of / is 
2ko/v = 2 and the initial value for h is 0.7. These agree with the initial values of the trajectory discussed in section 
II. 

The result is shown in Figs. [I] (continuous curves) and [i] We see that for small t the couplings / and h behave 
as in the model of section II, with / and h both increasing. At some point, however, A2 becomes sizable and then 
the last term in the right hand side of the second equation in |7]) pulls h towards zero. The trajectory is therefore 
characterized by a crossover from the IR regime where the fermionic interactions are mainly of Yukawa type, with 
IR free four-fermion interactions, and the UV regime dominated by the contact interactions, with UV free Yukawa 
coupling. This is similar to the behavior discussed in |llj . except that here we do not bosonize the contact interactions: 
we see that in the presence of the Goldstone bosons the crossover is automatic. 

IV. EXPERIMENTAL CONSTRAINTS 

In order to select a realistic trajectory among the various possibilities one would need to compare with experimental 
results. The model is consistent with EW precision measurements, as discussed in |12j . where the values of the EW 
parameters S and T are computed for the gauged NLerM; the addition of the fermion contact interactions does not 
modify this result. 




0.0 0.5 1.0 1.5 2.0 

t 

FIG. 3: A zoom at low-energy of Fig. [5] for the values (continuous curves) of the absolute values of the four-fermion operator 
coefficients A;. The dashed curve is the corresponding size for the contact interaction for the bound at A = 9.5 TeV. The dot 
at t = 1.55 on the dashed curve represents the bound at the reference value of the quark-level effective energy of 1.17 TeV. The 
values of the Xi must lie below the dashed curve for the experimental bound to be satisfied. 

Unfortunately the current bounds on contact interactions have been published only for the case in which a single 
operator, namely that proportional to our A3 is present [13] . This is rather unrealistic, given the RG mixing, but no 
need for a more detailed study was felt necessary. 

In this section we aim to show how the experimental bounds are in principle already able to tell us something about 
the size of the new interactions. Here we take, conservatively but rather unrealistically, the current experimental 
bound on A3 and enforce it on all coefficients as if they were contributing in the same manner to the partonic cross 
sections. In addition we assume that all three generations contribute with identical coefficients A^. A full Montecarlo 
simulation, and more comprehensive bounds, will be presented in another publication [14] . 
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The experimental bound is a lower bound of the so-called contact interaction scale A which is related to the 
coefficient A3 by the identity 

Mk) = -k\ (9) 

A bound on the value of A translates into a curve in the k, A3 plane because of the energy dependence in eq. As 
a crude estimate, we decide to impose the experimental bounds at the quark-level effective energy scale of the LHC, 
namely k e ff — \/s/(2 x 3) = 1.17 TeV, where the factor of 2 comes from the sharing of the energy with the gluons 
and the factor of 3 from the assumed equal energy partition among the three valence quarks. As mentioned before, 
for simplicity we will impose the same constraint to the absolute value of all four coefficients A^. 

As an example of the kind of constrains we can thus obtain, we take the most recent published bound [15] of A = 9.5 
TeV for an integrated luminosity of 36 pb _1 and values of A,; on a RG trajectory leading to the FP discussed in the 
previous section, namely fplc. Of course, the FPs with fewer relevant deformations lead to more predictions and 
therefore will be easier to disprove. As we can see from Fig. [3j the bound is satisfied. A future bound of A = 30 TeV 
is expected as the integrated luminosity of the LHC come close to 100 fb _1 |16j . 

In comparing with the experimental data, one may worry whether the power-law running of h and v may imply 
potentially important mass corrections to the cross sections. We compared the RG running of the masses for the 
solution above and verified that it does not deviate from the logarithmic one of dimensional regularization for more 
than 10%, at least below the scale so far explored, namely 600 GeV. 
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